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ABSTRACT 

We model the spectra and light curves of circumbinary accretion disks during 
the time after the central black holes merge. The most immediate effect of this 
merger is the deposition of energy in the disk due to the gravitational wave 
energy and linear momentum flux released at merger. This has the effect of 
perturbing the circular orbits of gas in the disk, which then intersect and radiate 
the dissipated energy. Because the disk is expected to be very optically thick, the 
radiation emerges predominantly in the infrared, and lasts for tens of thousands 
of years when the total black hole mass is M ~ IO^Mq. On the basis of a simple 
cosmological merger model in which a typical supermassive black hole undergoes 
a few major mergers during its lifetime, we predict that ~ 10^ — 10^ of these 
afterglow sources should be observable today. We also discuss the possibility of 
identifying them with existing multi- wavelength surveys such as SWIRE/XMM- 
LSS/XBootes and COSMOS. 

Subject headings: black hole physics - galaxies: nuclei 



1. INTRODUCTION 



Recent results from numerical relativity suggest that the merger of two rapidly spinning 
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). This recoil, or 



"kick," may lead to a number of interesting astronomically observable signatures, many of 
them much longer lived than the gravitational wave signal itself. Some of these signatures 
are indirect, and depend on galactic dynamics effects involving the merged black hole. For 
example, low-density cores in the central regions of galaxies may indicate a kicked black hole 



that ejected stars as it relaxed back to the center via dynamical friction (IMerritt et al.ll2004 
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Gualandris fc MerrittI 120071 ) . Alternatively, if the merged BH is ejected from the galaxy, it 



may also leave behind a bulge without chan ging its mass or velocity dispersion, thus breaking 
the M-o relation for a number of systems (jVolonterill2007l ). The overall population fraction 
of central black holes may g ive a way of indirectly measuring the distribution of BH recoil 
velocities (jSchnittmanl 120071 ) . 

Another class of recoil observations depends on the direct electromagnetic (EM) signa- 
ture from gas accreting onto th e kicked BH . The quasar that results may be spatially dis- 
placed from the galactic cen ter (Loebl 2007 ) , and it s line spectrum may be Dopple r shifted 



relative to the host galaxy (IBonning et al.l 120071 ). iMilosavljevic fc Phinneyl (120051 ) discuss 



the EM signal that would appear following a BH merger within a circumbinary gas disk, in 
which a central gap of gas was cleared out by the inspiraling BHs, but then is refilled on an 
inflow time after the merger. This scenario, independent of any recoil kick, could provide 
a soft X-ray counterpart to a LISA event on a timescale of a few months to years after the 
gravitational wave (GW) signal. More recently. iLippai et al.l (120081 ) investigated the emission 
from s hocks in a disrupted disk in a similar con figuration to th a t of iMilosavlievic fc Phinney 
(120051 ). but including a kick for the final BH. iKocsis fc Loebl (12008! ) estimated the (small) 
luminosity generated by viscous dissipation associated with the passage of the gravitational 
wave pulse itself through a surrounding accretion disk. Within this class of direct EM emis- 
sion models, new speculation predicts the possibility of t rails of gamma-ray ern ission from 
dark matter annihilation in the wake of the recoiling BH (IMohayaee et al.ll2008l ). 

In this paper we also consider a circumbinary disk around an inspiraling binary BH 
system, but focus on the long-term afterglow emission from the perturbed disk. At the 
time of merger, the gas in the surrounding disk is instantaneously placed on eccentric orbits 
due both to the change in gravitational potential from GW energy loss and the recoil of 
the central BH. These perturbed orbits relax back to circular trajectories while conserving 
angular momentum but dissipating energy over a few orbital periods. Because these disks are 
very optically thick, the dissipated energy will be radiated with a thermal spectrum, giving a 
warm yet relatively short-lived signal from the inner edge of the disk, followed by prolonged 
infrared (IR ) emission from the out er regions of the di s k. Th is signature differs from that 
predicted by lLippai et al.l (120081 ) and lShields fc Bonningl (120081 ) . who argued that there would 
be prompt emission in the ultraviolet (UV) and X-ray bands, based on calculations that 
do not include the disk opacity . O ur work also contras ts with previous efforts in that both 



Milosavljevic fc Phinneyl (120051 ) and ILippai et al.l (120081 ) focused on EM counterparts to LISA 



sources, and thus limited the total mass in the merger to < lO^M© (IKocsis et al.ll2007l ). We 
relax this limit because observations sugg est that the SMBH rnass function dN /d\og, M peaks 
closer to ~ 10^ Mq for redshifts z > 0.5 (IMarconi et al.l 120041 : iMerloni 1120041 ). However, we 
also consider the lower mass range so that we can estimate the character of the afterglows 
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produced by the smaller mass mergers that may produce gravitational waves detectable by 
LISA. 

For black hole masses of ~ IO^Mq, the IR afterglow signature could last for hundreds 
of thousands of years and provide evidence of SMBH mergers even without LISA signals. 
For reasonable pre-merger accretion rates, the central gap cleared out by the inspiraling 
binary does not close for ~ 10^ years, giving a unique signature with infrared luminosity 
comparable to that of a quasar, and virtually no UV or X-ray emission. Similarly, unlike 
classical AGN, we expect no compac t radio je t s or p hotoionized emission lines. Using the 



anti- hierarchical BH ma ss function of iMerloni I (120041 ) and an estimate of the SMBH merger 



rate (ISesana et al.ll2004l ). we expect a total merger rate of ~ 0.1 per year out to z ~ 6 for 
total mass M > IO^Mq. Due to the long lifetimes of such afterglows, as many as ~ 10^ — 10^ 
may be observable in the entire sky at any one time. Within the fields and flux sensitivities 
of existing multi-wavelength surveys such as SWIRE and COSMOS, we expect that ~ 1 — 10 
may be discernable today, but any candidates would require follow-up optical spectroscopy 
for positive identification. 

This paper is organized as follows: in Section [2] we describe a simple model for a binary 
merger including GW energy and momentum losses and the subsequent behavior of the 
perturbed circumbinary disk, with predictions of light curves and spectra. In Section [3] we 
combine these results with a cosmological merger model and calculate expected merger rates 
and distribution of afterglow spectra. In Section H] we discuss the feasibility of observing 
such systems, and present our conclusions in Section O 



2. SPECTRAL SIGNATURE OF A SINGLE MERGER 

2.1. Circumbinary disk model 

We begin with a circumbinary disk described by an inner radius -Rin, total BH mass 
M = Ml + M2, mass ratio q = M1/M2 < 1, accretion rate M, and an efficiency rj such that, 
if the disk extended all the way in to a single central black hole, the luminosity would be 
L = TjMc^. For this efficiency, we can define a normalized accretion rate fn = M/M-^dd, 
where LEdd = V^EddC^ is the Eddington luminosity. Outside of -Rin, we describe the disk by 
a steady-state alpha mode l, with nominal efficiency 77 = 0.1 and viscosity parameter a = 0.1 



(IShakura &: Sunyaevlll973l ). The sensitivity of our results to these assumptions is discussed 



in Sectional 

For a pair of BHs on a circular orbit with binary separation a, the inner edge of the 
circumbinary disk is located at Ri^ = Xa, where A(g) is a function only of the mass ratio and 
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exhib its a relatively small range around A ^ 1.6 — 1.8 for q > 0.01 (lArtymowicz fc Lubow 
19941 ). Interior to this point, gas either accretes directly onto one of the two BHs or, 
more likely, is ejected from the s ystem, giving a gap o f very low density inside of -Rjn 
( iMacFadyen fc Milosavljevidl2008l : iBogdanovic et al.ll2008l ). When the binary separation is 
sufficiently small, its evolution becomes dominated by gravitati onal wave los ses and inspirals 
on a timescale given by the leading-order quadrupole formula (jPeterslll964l ): 



t 
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where ^ = M1M2/M is the reduced mass of the binary. 



As the binary orbit shrinks, the inner edge of the disk follows closely behind, maintaining 
the relation i?in = Xa as long as the inflow time is less than th e inspiral time. For a standard 
alpha disk, the gas inflow time is given by [e.g. iKrolikI (119991 )] 



^mfiow(-R) 



orb 



^3/2 



act 



a 



(2) 



where VLorh is the orbital angular velocity, Cg is the sound speed, and h is the disk thickness. 
In the inner regions of the disk, where the pr essure is radiation-domina t ed and the opacity 
is do minated by electron scattering, we have (jShakura &: SunyaevI Il973l : iNovikov fc Thorne 
I973I ) 

tinflow(i?rad) ^ 5 X 10*^ aZ\rj\mZlMsxy^ yr, (3) 
while in the outer, gas-pressure and scattering-dominated region, we have 



yr. 



(4) 



Here we have employed scaled parameters a = 0.1q;_i, t] = 0.1?7_i, m = 0.1m_i, M = 
IO^MsMq, and x = 10^ x ^iRc^ /GM) is the radius in geometric units. The point where 
the disk transitions from radiation- to gas-pressure dominated can be estimated by setting 

^inflow (-^rad) ^inflow (-Rgas) ; OT 



Xtr 



4 X 10^ ^,^10/105 



Scaled in terms of xs, the inspiral time is 



tinsp(i?in/A) ^ 1 X 10*^ Ms 



■A *{q)xi yr. 



(5) 



(6) 



When the inspiral time becomes shorter than the gas inflow time at the inner edge, the 
binary separation decreases faster than the gas can move in, effectively decoupling the BHs 
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from the disk, and the merger occurs soon after (jMilosavljevic fc Phinneyl l2005l ). For our 
nominal model parameters M = 10^ Mq, g = 1, 77 = 0.1, and rh = 0.1, the central gap in the 
disk is quite large: -Rin ~ lO^M. However, we also find that for some system parameters, 
namely when q is small (thus a long inspiral time) and the accretion rate is large (thus a 
short inflow time), the inner edge of the disk is able to keep up with the gravitational inspiral 
up to the final plunge and no appreciable gap is formed. Furthermore, we find that whenever 
a gap is formed, the inner edge is almost always within the gas-pressure dominated region 
of the disk, so we can solve for R^^ with equations (jH E]) : 

5/13 



L(l + g)2 



(7) 



(Note that in iMilosavljevic fc Phinneyl (120051 ) . the inflow time is taken to scale as R^, and 
they consequently find smaller values of Rm, typically within the radiation-dominated region 
of the disk.) 



2.2. Post-merger dynamics of a perturbed disk 



At the time of merger, M is reduced nearly instantaneously by a fraction ecw due to 
GW energy flux, with ecw typically a few percent. This mass loss is scaled to the symmetric 
m ass ratio u = q/(l + q )^i normalized to nume rical simulations of non-spinning BH mergers 



(IGonzalez et al.ll2007al : iHerrmann et al.ll2007al l: 

ecw ~ 0.08 z/ + 0.32 1/2. 



(8) 



This mass loss could be as high as ecw ~ 0.1 for aligned, rapidly spinning BHs (IMarronetti et al 
20071 : iDain et al.ll2008l ). but averaging over an ensemble of uniformly distributed spin orien- 
tations, we think the non-spinning estimate is reasonable. 

For generic systems without special symmetries, the resulting black hole also receives a 
linear momentum recoil with velocity Vkkk- This recoil leaves the outer regions of the disk 
(where Vorh ^ M^ick) unbound and significantly disrupts the inner regions. For the bound 
regions (Krb ^ ^kick), an annulus of mass dM{R) = 27rRdRT,{R) on an originally circular 
orbit at radius R receives an instantaneous boost in energy relative to the merged black hole 
of 

,^ , GdM(R)M 1 ^ „ 



R 

In the unbound outer regions, the change in energy is simply the original binding energy 
dEo: 

GdM{R)M 1,^^,^,^. 



C^-E'kick(-Runbound) — —dE^ 



2R 



-dM{R) K^b- 



(10) 
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In the inner-most regions, where the mass-loss contribution to the energy gain is sig- 
nificant, it is possible that not all this added energy will be available for dissipation. This 
is because (as pointed out to us by Cole Miller), in the perfect (i.e. dissipationless) fluid 
limit, fluid elements retain their initial angular momentum. With the new, smaller central 
mass, that angular momentum is too large for any bound orbit at their initial radius. If fluid 
dynamics enforces circular orbits in the long-run, individual fluid elements expend most of 
the energy they gained by the mass-loss on moving out to their new orbital radius, leaving 
only ~ {GM/R)eQ^r available for dissipation. However, the actual mechanics of this response 
may well involve numerous shocks in which considerable angular momentum mixing between 
different fluid streams can occur. The effectiveness of this mixing will depend on numerous 
considerations including the radial surface mass density profile, the rapidity of cooling, and 
the role of radiation pressure (see below for further estimates relevant to the latter two), and 
its evaluation is therefore best left to future work. Given this uncertainty, we have opted for 
the more optimistic scaling (heating from mass-loss ~ ecw times the binding energy) in the 
results shown here. However, we have checked that whether one adopts a mass-loss heating 
scaling ~ ecw or ~ makes very little difference to our conclusions because it alters the 
afterglow at an interesting level only during the first 10^-10"^ yr after merger, while nearly 
all the detectable objects will be seen at much later times. 

For the systems we have considered, the vast majority of the afterglow energy is released 
in the regions of the disk that are dominated by gas pressure (typically R > lO^M ), where 



the surface density is given by fIShakura fc SunyaevI 1 19731 : iNovikov fc Thornd Il973l ) 



S(/2) ^ 1 X 10'^aZt^%zTmJ^M'J''xf^'' gm/cm^ (11) 

and the gas density is 

p{R) ^ 1 X 10-'^aZl^'%ZlVJ^M-'/'\f'/'' gm/cm3. (12) 

From here on, we use a fiducial radial distance of 10"^ gravitational radii, corresponding to 
the marginally bound region of the disk for the largest expected kicks. Strictly speaking, 
equations flll^ [T^ assume that the opacity is dominated by electron scattering, but they are 
actually quite similar to the expressions corresponding to free-free/bound-free scattering. In 
practice, both opacity effects tend to be important (and also dust and molecular opacities at 
lower temperatures in the outer-most regions of the disk) and we find our pre-merger disks 
are extremely optically thick out to radii of at least x = 10^. Since the total energy released 
at each point in the disk is proportional to the surface density (eqn. [H]), only optically 
thick disks will actually have enough mass to produce a significant amount of luminosity. 
Moreover, because S(-R) depends on rh only to the 3/5 power and on M only to the 1/5 
power, our conclusion that these disks have large optical depths should be valid over a wide 
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range of possible conditions. As we will show below, even the abrupt heating that takes 
place after the black hole merger does not alter this conclusion. 

The change in energy given in equation leads to a perturbed disk with eccentric, 
intersecting orbits. Conserving angular momentum in each annulus, the gas relaxes back 
to a collection of circular orbits on a timescale of order theat as shocks due to intersecting 
orbits dissipate energy in the disk. We define the heating time as the time it takes for nearby 
eccentric orbits to cross. In the epicyclic approximation for small eccentricity e, the radial 
coordinate i? of a perturbed orbit is 

R{t) = Roil + ecos-^), (13) 

where Rq is the radius of the guiding center orbit and = ^orb^ + "^o is the epicyclic phase. 
Density caustics form when orbits with different Rq overlap at the same R: 

dR 2{GMf'H, . M M n 

-mr = o ^7^(esmV') + (l+ecos7/;) =0. (14) 

In the limit of small e, we get a heating time of 

, ^ 2 i^O^ _ 1 torb , . 

In the simplest case where the kick velocity is zero, the instantaneous GW mass-loss will 
excite the entire disk with a constant eccentricity Cmass-ioss ~ cgw; independent of radius. 
In the opposite regime where the mass-loss contribution to the eccentricity is much smaller 
than that of the recoil, the post-kick magnitude of the total specific angular momentum ^ of 
a single fluid element relative to the merged black hole will be 

e = Rl{V^,^ + 2/3KrbVkick + /5VkL + ^VkL), (16) 

where /3Vkick is the component of the kick parallel to the fluid velocity vector, -zVkkk is the 
component of the kick perpendicular to the orbital plane. The semi-major axis of such a 
fluid element can be determined from its specific energy e\ 

Combining equations fll6|17l) . we find the perturbed fluid element will have a post-kick 
eccentricity of 

(1 + 3/52 _ ^ ^2/5(1 + + z^) + + z") 



^kick 



1/2 

I^kick 



GMa / K 



orb 
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This is an exact result for bound, Keplerian orbits. For V^ick ^ K^rb, we find the typical 
eccentricity scales like e ~ Vkick/Krb for planar kicks with z = and e ~ (Vkkk/Krb)^ 
for kicks directed out of the orbital plane with z = 1 and (3 = 0. While the final kick 
directi on is not independent of its magnitude, the relationship is not well-known at this 
point (ITichy fc Marronettil 120071 ). so for these simple estimates, we typically average over 



angles to get 



1/3. 



In practice, when estimating the heating rate in equation f[T^ . we use a combination of 
the mass-loss and kick eccentricities: 



etot 



=GW 



+ {e 



2 ) 

kick/' 



(19) 



where {e^ick) the average of equation (fTSll over an isotropic distribution of kick directions. 
We further impose the constraint that etot < 1 everywhere in the disk, with etot = 1 indicating 
the unbound regions. 

To calculate the light curves produced by the merger, we must first estimate the tem- 
perature, density, and optical depth of the perturbed disk. During the heating phase, shocks 
with characteristic speed Khock will be driven through the gas at each annulus in the disk. 
In much of the bou nd region, the shock speed will be roughly limited by the kick speed, a 
result confirmed by IShields fc BonningI (120081 ). who calculated the relative velocity of col- 
liding geodesic particles on perturbed orbits, finding Vshock ~ O.SVkkk at early times and 
^hock ~ O.Q Vkink at later times. In th e innermost regions of the disk, the GW mass loss 
(neglected by [Shields &: BonningI (120081 )) will dominate the orbital dynamics, giving 



^hock, bound < max [(2eGw) Krb, ^kick] • 

In the unbound regions, the shock speed is hmited by the orbital velocity, so 

^hock.unbound ~ min [V"orb; H;ick] 



(20) 



(21) 



everywhere in those regions. It should be noted that even if the s hock speeds i n the b ound 
regions are much smaller than the kick speed, as suggested by iLippai et al.l (120081 ) . the 
total kick energy must still be dissipated eventually, and only the initial temperatures will 
be different. However, as we will see below, the post-shock disk should be quite optically 
thick to thermalization, giving robust spectral and light curve predictions independent of 
the specific heating details. 
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2.3. Dissipation of energy in the disk 

Assuming the pre-merger disk is relatively cold, the initial shocks heat the gas to a 
temperature 

^2 



-t shock — shock -"^ 



1.4 X 10^1^000 K, (22) 



where /x ~ 0.6mp is the mean molecular mass of the gas and Viqqq is the shock velocity in 



thousands of km/s (IMcKee fc Hollenbachlll980l ). This energy in turn will be transformed to 



radiation on the bremsstrahlung radiation timescale: 



, _ ^kTsl^ock ~17T/ n.^/^%2/5rh-2/5 /I #7/10 33/20 .^oX 

trad — — — : ^ ~ -L Kiooo Ct-l V~l^~l ^^-^S ^4 l^'^J 



1.4 X 10~^'^ip/m,)T:^:^^ 



where we have used the free- free volume emissivity given in IShapiro fc TeukolskyI (119831 ). 
This timescale is much shorter than any dynamical time in the system, so we assume all the 
shock energy is converted instantly into radiation. 

At the high temperatures prevailing immediately post-shock, the disk may no longer be 
optically thick to free-free absorption, but will still be highly opaque to electron scattering. 
Where, as here, Tes > Tff, the optical depth to thermalization is given by rtherm = (TcsTfr)^^^, 
where Tes and tq are the electron scattering and free-free optical depths, respectively: 

= 0.35 E, (24a) 
%S ^ 2.3 X 10^^ pT-^/^j^, (24b) 



Tthcrm - 0.5a:f/^^:fmi/,^M-^/^^,-^^/^V-^/^ (24c) 

Here Rg is the mean free-free opacity calculated by averaging over a bremsstrahlung emission 
spectrum at temperature T = 10^ Ty K. Thus, for typical parameters, we expect the disk 
immediately post-shock to be marginally optically thick in terms of thermalization. 

When the parameters are such that the initial burst of free-free radiation is well- 
thermalized, the spectrum quickly evolves through absorption and reradiation to Planck 
form. As it does so, the matter and the radiation reach a state of thermodynamic equilib- 
rium at a temperature defined by matching the thermal energy density created by the shock 
to the total heat content of gas and photons: 



(25) 
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For any shocks with \4hock ^ 100 km/s, the thermodynamic equihbrium pressure will be 
completely dominated by the radiation. The equilibrium temperature is then 

T,,_ ^8x10^ l/,Voo«:r^:y^°--/i °M3-^/^\l^^/«° K. (26) 

The fact that the thermodynamic equilibrium is strongly dominated by radiation pres- 
sure has the consequence that only a small portion of the initially-radiated free-free emission 
need be immediately thermalized in order for the gas quickly to reach thermodynamic equi- 
librium. Because free-free opacity at frequencies v well below kT/h is oc even when the 
blackbody peak is at first not thermalized, sufficiently low frequencies are. Thus, for those 
regions with Ttherm < 1, we can define a cutoff frequency t'therm below which the spectrum 
will be thermalized: 

{ ^ 1.3 aZf''\-_Tm'llM-'''\l'''T-''\ (27) 

V "^-t shock / 

If Ttherm < 1, the thermafized band will always be in the Rayleigh- Jeans portion of the Planck 
spectrum. 

The total energy content in this low-frequency thermalized portion of the spectrum is 
given by 

ttherm ~ ^ ( , r^""^ ) Ct^shock- (28) 

When this is comparable to the energy in the post-shock gas, 3/2{p/ fi)kTshock, the tem- 
perature of the matter must be reduced by a factor of order unity in order to supply the 
energy in the radiation, in turn raising Ttherm (fx T"'^/^) and therefore f/therm- Further 
depression of the gas temperature follows, leading to a still wider bandwidth of thermaliza- 
tion. Rapid approach to thermodynamic equilibrium is thus the end-result whenever the 
initial energy in thermalized photons is comparable to the gas energy. This condition of 
ttherm ^ 3/2 (p//i) /cTghock wiU be satisfied for all shock temperatures with 

T.hock < 2 X loVlVTrh'/^Mi^W^' K. (29) 

For our fiducial parameters, the corresponding shock speed is so high — nearly 40, 000 km/s, 
roughly ten times the largest possible kick speed — that highly unusual conditions are required 
for this condition to be violated. Furthermore, since the shock speed is limited by the 
orbital speed, and V"orb ~ 3000^4 km/s, only systems with extremely low accretion rates 
{rh^Ji ^ 10"^) have any chance of avoiding rapid thermalization. But as we mentioned 
above, since the luminosity is proportional to the disk surface density (and S ~ rfi^^^), we 
expect to observe only those systems with moderate-to-large rh in the first place! 
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Even if free-free absorption were not enough to thermalize the radiation, Compton 
scattering could also do the job on a timescale roughly comparable to trad- This is because 
the characteristic Compton-?/ parameter at the post-shock temperature is 

y = ^-^rl ^ 1 X io« r,^ooo«:f (30) 

In this case, the thermal energy of the gas is drained as the electrons scatter into a Wien 
spectrum the large number of low-energy photons they have created by free-free emission. 
The timescale for this process is almost as fast as the free-free radiation cooling time itself be- 
cause each low-energy photon receives on average an additional energy ~ kT^ by subsequent 
inverse Compton scattering. 



2.4. Light curves and spectra 



Thus, by any of several processes, the very large column densities found in these disks en- 
sure that thermalization is inevitable. After the disk is shock-heated and reaches a radiation- 
dominated thermal state, it will expand adiabatically to a post-shock scale height determined 
by hydrostatic equilibrium: 



^shock ~ 



'^s orb ~ 



with the pressure and density given by p = (i£'kick/(27r/ishock -Rc^-R) and p = S//ishock- Since 
the initial orbital energy of the annulus is (IEq = (27TRdRT,){—GM/R), we can write the 
post-shock scale height as 



h 



shock 



dEn 



1/2 



R. 



(32) 



The internal energy of the perturbed disk is then radiated over a cooling time at each radius, 
which for optically thick, radiation-pressure dominated disks is given by 



tcool -Tes- ~0.33E(i?) 

c 



dE, 



kick 



dEn 



1/2 



R 



(33) 



The evolution of internal energy dEghock due to intersecting orbits of the perturbed disk 
is thus governed by 



d 

"77*^^ shock 

dt 

which gives a radiated luminosity of 

4^ _ dE^ick{R) 

"-t^shockl-K; fcj — — 7— X 

thoatl-Kj 



dEv, 



kick 



^hcat 



dE. 



shock 



^cool 



/tcool l^je"*/* 

cool 



t < ^heat 
t > ^heat 



(34) 



(35) 
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Thus, the luminosity grows hnearly at early times, which is reasonable for crossing eccentric 
orbits, and decays exponentially at late times, appropriate for a diffusion process. 

On top of this luminosity from the internal heating of the perturbed disk, we also assume 
that outside of Ri^, the disk continues to behave as in a steady-state accretion disk, removing 
energy and angular momentum from the gas, giving the classical accretion luminosity at each 
point in the disk: 

dLacc(i?) = 7;—F^ 5-"^- (36) 

While the shock-heating and subsequent expansion of the perturbed disk may change the 
local accretion rate, we generally find the cooling time to be much shorter than the inflow 
time, so we continue to assume an overall steady-state disk model with constant M at all 
radii. We anticipate that more detailed numerical simulations should be able to model the 
dynamic behavior of these perturbed disks in the near future. The local emission spectrum 
is taken to be a thermal blackbody with 

dLtot{R) = c/i^shock + dL,,, = AttR dR aT\R), (37) 

with a the Stefan-Boltzmann constant. After the merger, the inner edge of the disk con- 
tinues to migrate inwards on an inflow timescale, eventually closing the gap and forming a 
"standard" accretion flow in the inner regions of the disk. 

Figure [U shows light curves for several events with M = IO^Mq, t] = 0.1, and a = 0.1, 
but varying q, ^kick? and rh. In all cases, we see an initial linear rise in luminosity as the gas 
near i?in is first excited and then radiates its heat on a timescale of a few years. The peak 
luminosity is typically reached after ~ 10'^ yr, and the afterglow remains roughly this bright 
for ~ 10^ yr. In the solid, dashed, and dot-dashed curves, the binary mass ratio is unity, 
giving an inner edge of -Rin ~ lO'^M. In these cases, the inner gap is not closed until nearly 
10^ years after merger. At that point, the luminosity increases sharply. The reason for this 
abrupt behavior is that the inner edge of the disk moves in very slowly at large R, then 
accelerates its inward motion as the gap closes, leading to a nearly instantaneous increase 
in total luminosity. The solid and dashed curves have m = 0.1, while the dot-dashed curve 
has rh = 0.2. For a higher accretion rate, -Rin moves in slightly, and from equation fill I) , we 
see that overall mass-and thus luminosity-of the perturbed disk also increase with m. 

The solid and dot-dashed curves in Figure [1] correspond to a recoil of Vkick = 1000 km/s, 
and the dashed and dotted curves have Vkick = 300 km/s. From equation we see that for 
the inner regions of the disk, where GM ecw/ R 3> V^^^k; perturbations due to mass loss 
are greater than those due to the kick, while in the outer regions, the kick term dominates. 
The dotted curve has q = 0.1, which causes the inner edge of the disk to be closer in at 
the time of decoupling, so the delay before gap filling is also shorter. From equation ([8]), we 
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expect a smaller GW mass loss, and so the initial heating rate is smaller. As can be seen 
from equation fl35|) . the inner disk lights up on roughly the same timescale as the g = 1 cases 
despite having a smaller value of -Rin- Note that all cases eventually settle down to nearly 
the same luminosity after the gap closes because the luminosity in the normal AGN phase 
is determined only by rh and M. 

In Figure [2] we show a time sequence of spectra from a merger with M = 10^ Mq, 
m = 0.1, g = 1, and l^ick = 1000 km/s (the solid curve in Fig. [T]). Over the first hundred 
years or so, the disk brightens as the perturbed orbits near the inner edge begin to intersect 
and dissipate energy. Then, as the region of maximum dissipation propagates outward, the 
disk slowly dims and reddens, its radiation moving from the optical/UV to IR over thousands 
of years. As in Figure [H at around 10^ years after the merger, the inner gap closes, forming 
a "typical" AGN with thermal emission peaking in the UV band. 

At even later times, the situation becomes highly uncertain, as the system's behavior 
depends largely on the existence and state of gas in the outermost regions of the disk {R > 
lO^M). Even within this r adius, it is not clear whether the gas can be described by a 



simple thin disk model as in IShakura fc SunyaevI (119731 ). Beyond R > lO^M, the disk may 
become gravitationally unstable, but there appears to be significant observational evidence 
that complete fragmentation and collapse is avoided by some auxiliary heating mechanism 



ipj. 

Lodatol (120081 ) and references therein]. In our simple model, if there is in fact appreciable 
gas and thus emission beyond R ~ lO^M, it should appear as a steadily reddening peak in 
the far infrared, as might be extrapolated from the curves in Figure [2J However, even if 
there is a significant amount of gas in these outermost regions, it will almost certainly be 
unbound after any appreciable kick (orbital speeds at i? ~ lO^M are a few hundred km/s), 
and thus inherently limited in its total available energy content. 

To estimate the light curves and spectra we might expect to see from LISA counterparts 
with M ~ lO^M©, recall that the total energy released in the afterglow scales like Si? ~ 
^Yie typical timescales vary as M. Thus the peak luminosity should be a factor of 
~ 10~^ smaller than those plotted in Figure [1] and evolve on a timescale 100 times shorter. 
The surface brightness, however, which scales like S, has only weak dependence on the total 
mass, so we expect LISA afterglows to have similar spectra, peaking in the near- to mid-IR 
over a period of a few years after the merger. 
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COSMOLOGICAL MERGER HISTORY 



We adopt a standard ACDM cosmology w ith = 0.74, Qm = 0.26, and normalized 
Hubble constant /i = 0.73 fISpergel et al.l 120071 ). At each red shift, we approx imate the BH 
mass function by the anti-hierarchical distributions given by iMerloni I (120041 ). extrapolated 
to redshift z = 6. We assume that each BH experiences an average of iVmm major mergers 
between z = 6 and the present an d that over this rang e, the merger rate per comoving volume 
is weakly dependent on redshift (ISesana et al.ll2004l ). The top panel of Figure [3] shows the 
expected number of mergers with M > IO^Mq throughout the universe per observer year 
per unit redshift fo r A^mm = 3 (in rough agreement with the lower-right panel of Fig. 1 in 
Sesana et al.l (120041 )). All rates quoted below should simply scale linearly with this parameter 
A^mm- Note that these merger rates are significantly lower than those typically quoted for 
LISA sources, wh ich include a large number of smaller black holes (M < lO^M fr^) at high 
redshift [z > 10) jMenou et al.ll200ll : ISesana et al]l2004l : iRhook fc Whvithelboosh . 



While the mass functions given by iMerloni I (120041 ) likely underestimate the number of 
low-mass AGNs due to observational selection effects, these systems are also much less likely 
to produce long-lived, bright afterglo w signals. There fore, in estimating number counts, it is 
reasonable to use the observations of iMerloni I (120041 ) to normalize the high-mass portion of 
the distribution function, even though they may be interprete d in te rms of anti- hierarchical 
evolution, and use the hierarchical simulations of ISesana et al.l (120041 ) to estimate the merger 



rates of these high-mass systems as a function of redshift. As mentioned above in Section [2^^ 
even if the number of low-mass (M < lO^M©) binaries is greater by an order of magnitude, 
their luminosities and lifetimes will also be much smaller than those systems with mass 
around IO^Mq, so the observable IR luminosity function should be rather insensitive to 
mergers involving low-mass black holes. 

For a given merger, we determine t he masse s by se lecting two black holes randomly 
from the distribution $(M, ^) giv en in iMerloni I (l2[)04r) . The kick magnitude is deter- 



m ined from the formulae g i ven by Baker et al. ( 20081 ). which agree with earlier estimates 



of iSchnittman &: Buonannoi (120071 ) for small to moderate kicks, but predict a larger number 
of extreme kicks above 1000 km/s. We assume a uniform distribution of BH spin orientations 
and set a/M = 0.9 for all systems. Recoil velocities generally scale linearly with the spin 
parameter, so should not be t oo sensitive to small unc ertainties in this typical value. With 
regard to the spin orientation, iBogdanovic et al.i (120071 ) recently argued that BHs embedded 
in a circumbinary accretion disk should rapidly align their spins with the overall angular 
momentum of the disk, which would in turn lead to much smaller recoil values. On the other 
hand, if the disk forms outside of the binary, we expect very little gas to actually accrete 
directly onto either black hole, and the BHs could very well maintain their original random 
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orientations fISchnittman II2004I ). 



While the vast majority of isolated SMBHs are not active at any given time, the accre- 
tion rates for mer ging BHs may oii average be much higher due to the inflow of gas afte r 



galactic mergers (IHaehnelt fc Reed Il993l : iMihos fc HernquistI 1 1994 IVolonteri et al.l l2003l ). 



Furthermore, many models rely on the dynamical friction of gas disks to bring the BHs 



(Escala et al. 


2005; 


Dotti et al. 


2007) 



ing system has a significant circumbinary disk, with accretion rates uniformly distributed in 
logarithmic space between m = 0.01 and 1. 

For each merger, we define the afterglow phase as the time during which the total 
luminosity in equation 0371] is dominated by the shock heating due to the perturbed orbits 
and where the inner gap has not yet closed. In the bottom panel of Figure [3] we show 
the expected number of objects in the afterglow phase observable at the present time as a 
function of redshift. Note that while the merger rate (top panel of Fig. [3]) decreases with 
z > 1, the redshifted time makes those events appear to last longer to us, and thus there 
should be more in the afterglow phase at any given time. Additionally, due to the anti- 
hierarchical growth of SMBHs, mergers at higher redshifts typically involve larger masses 
and longer afterglows. 

In FigureHlwe show the luminosity function $giow(-^, z) (in units of number per comoving 
volume per log luminosity) for systems in the afterglow phase in a number of different redshift 
bins. The distribution is remarkably narrow and constant in time. There are a number of 
reasons for this behavior. First, the BH mass dist ribution we use is peaked around lO^M© 
and falls off sharply above lO^M© (IMerloni 1120041 ). In the downsizing paradigm of cosmic 
evolution, the number of smaller-mass BHs (< lO^M©) increases at low redshift, in turn 
giving a slightly larger number of low-luminosity afterglows, as shown by the black and blue 
curves in Figure HI However, these low-mass systems are also much shorter-lived in their 
afterglow phase, giving a strong selection effect against seeing them at any one time. The 
mass accretion rate rh is also limited by a selection effect: we consider only rh > 0.01 in our 
model because if rh is too low (thus giving a low-luminosity, low-mass disk), the BH binary 
does not evolve quickly enough via dynamic al friction with the disk to merge within a Hubble 
time (lEscala et al.l l2005l : iDotti et al.l 120071 ). If, however, rh is too large (high- luminosity), 
the inflow time at the inner edge of the disk [eqns. ([3l H])] will be small enough such that no 
appreciable gap is formed and the afterglow is not readily identifiable. Additionally, we see 
from equations and (ITTi) that the total energy released is a relatively weak function of 
rh, further narrowing the range of luminosities. Lastly, the mass ratio q is likely to be close 
to unity for a similar reason: if it is too small, the GW inspiral time is longer than the gas 
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inflow time and the gap is small and short-lived. This selection for large q also favors higher 
kick velocities, which fo r spinning BHs are maximized when q = 1 (jSchnittman fc Buonanno 



20071 : iBaker et all 120081 ). 



4. OBSERVATIONAL POTENTIAL 

Because any single merger event spends the vast majority of its lifetime in the late, 
IR-dominated phase corresponding to the outer disk relaxation, the average system seen 
today should be characterized by the double-hump spectra shown as the orange and yellow 
curves in Figure [H Using the cosmological model parameters from the previous section, 
in Figure [5] we show a randomly selected sample of rest-frame spectra within ~ 10^ years 
after merger. The solid curves correspond to those systems formally in the afterglow phase, 
when the luminosity is dominated by the dissipation from the kicked disk. For reference, 
the dashed curve is from a system that has already closed its central gap and radiates as 
a normal AGN in the inner regions. We find that, for the model parameters used above, 
the numbers of closed-gap and open-gap systems are roughly equal, with somewhat more 
closed-gap systems at low redshift, where the typical BH masses — and thus timescales — are 
smaller. 

Closed-gap systems are essentially identical to "normal" AGN, dominated by a thermal 
peak in the UV from the inner disk, along with a strong IR peak due to reprocessing of 
the UV by a surrounding dusty torus at large distance, which could easily be confused with 
the IR afterglow emission from the perturbed outer disk. Thus, for identification purposes, 
we will focus on kicked disks where the central gap has not yet closed, so the emission is 
almost entirely in the IR. Due to the lack of a central engine, these systems are not likely 
to produce significant UV/X-ray fiux or compact radio jets. Furthermore, we expect their 
time variability to be quite low, since most of the emission is coming from the outer regions 
of the disk, where the dynamical time is hundreds of years. 

Restricting ourselves to this sample of open-gap disks, in Figure M we show the source 
counts in the sky over a range of observed wavelengths, plotting contours of N{> S^, Xohs)- 
Also shown in Figure [6] are the 5a fiux limits from the wide-field, multi-band SWIRE and 
COSMOS surveys. SWIRE has a total coverage of nearly 50 deg^, so on the basis of Figure 
[5] we might expect to find ~ 10 afterglowing sources in the SWIRE field. The COSMOS 
survey covers a somewhat smaller area, 2 deg^, with sensitivity comparable to SWIRE at 
24 fim, but significantly better at 3.6 and 8 /im, so there may again be ~ 1 — 10 afterglow 
sources detectable in its field. The GOODS survey goes considerably fainter, to 40 /iJy at 
24 /im, but covers an area of just 150 arcmin^, giving an expected number of sources of only 
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~ 0.01 - 0.1. 

In order to distinguish these afterglow signals from other ultra-luminous IR sources, we 
require very low flux limits from optical, UV, and X-ray observations. As an initial estimate, 
we might require an X-ray flux no greater than 10% of the total IR flux, corresponding to 
roughly 10~^^ — 10"^'^ ergs/cm^/sec for the typical afterglowing source within z <3. This is 
well within the limits of the XMM-Newton Large Scale Structure and the Chandra XBootes 
surveys, which together cover roughly 18 deg^ of the SWIRE field. COSMOS has somewhat 
deeper coverage with XMM, as well as full HST coverage with ACS, so may be the preferable 
approach at this point. 

As we can see from Figure O in the mid-IR, the A^(> S) curves scale roughly as S~^, 
giving a larger number of high-fiux sources than a uniform distribution with no cosmic 
evolution (N ~ S~^^'^). Thus the number of potentially confusing sources should decrease 
with increasing flux, yet detection will require covering an observed solid angle that increases 
linearly with flux. In light of these trade-offs, we believe the best discrimination potential 
may ultimately come from an extremely wide-field survey, even if shallow. To estimate the 
distance and brightness of the high- flux end of the distribution, in Figure [7] we show source- 
count contours as a function of redshift, finding that the majority of bright sources are indeed 
nearby, with 2; ^ 1. 

Many IR/X-ray selected candidates are likely to be obscured AGN with a high Compton 
depth, which may be ruled out as afterglow sources by the detection of hard X-rays or 
compact radio jets. In addition to demonstrating a low level of X-ray flux, we may need 
optical spectroscopy to rule out the existence of line emission from UV excitation of the 
region around the AGN. This may require targeted spectroscopic observations, focusing on 
candidates selected with the other criteria described above. 



5. DISCUSSION 

We have proposed the existence of a new class of electromagnetic sources produced by 
the merger of two supermassive black holes and thus corresponding to some of the strongest 
gravitational wave signals in the observable universe. These objects are potentially observ- 
able today, long before any planned GW experiment might detect them. When surrounded 
by a circumbinary accretion disk, the two BHs can be driven towards merger via dynami- 
cal friction. Eventually the gravitational wave losses dominate the evolution of the binary, 
shrinking the orbit faster than the gas inflow time, at which point the disk decouples from the 
binary, leaving an open gap of extremely low-density gas and negligible direct accretion onto 
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the BHs. This decouphng is followed soon after by the merger of the two BHs, disrupting 
the gas disk by the energy and momentum losses in the gravitational waves. The perturbed, 
eccentric orbits then intersect, forming shocks and heating the disk, which then radiates 
its internal energy over a cooling time. The large optical depth in the disk ensures that 
however the heat is released initially, it will rapidly reach thermal equilibrium and should 
therefore emit primarily in the infrared over most of the afterglow's lifetime. Furthermore, 
only systems with high optical depth will also have sufficient mass in the disk to produce 
significant luminosities from internal shocks. Due in part to the large optical depths of the 
disks, these IR afterglow signatures could last hundreds of thousands of years for typical BH 
masses of M ~ 10^ — IO^Mq. At some point long after the merger, the inner edge of the 
disk will migrate in towards the central BH, forming a normal AGN and marking the end of 
the afterglow phase. 

Folding this perturbed disk model with a cosmological merger scenario, we predict the 
expected event rate (~ 0.1/yr) and total number counts (~ 10^) for afterglowing disks ob- 
servable in the universe today. Considering only existing survey data from Spitzer, XMM, 
and Chandra, roughly 1 — 10 open-gap afterglow systems may be detectable today. However, 
these sources would need to be discriminated from a much larger population of deeply ob- 
scured AGN, either with hard X-ray observations or optical spectra that might rule out the 
existence of narrow-line regions characteristic of normal AGN. With a large fraction of their 
luminosity coming from the outer regions of the disk, where the dynamical time is hundreds 
of years, we also expect to see relatively low levels of variability in their light curves. 

In addition to the observational challenges in successfully detecting and identifying these 
systems, there are also a number of theoretical uncertainties in the underlying afterglow 
model. Most critically, the properties of the gas disk (if it even exists!) at large radii will 
strongly affect the predicted light curves and spectra described in Section [231 At this point, 
we have relatively little understanding of AGN disks beyond R ~ lO^M, where much of the 
afterglow emission originates. Even if the disk reliably extends out to R > lO^M, the gas 
accretion rate as well as the surface density and temperature in these regions may very well 
not be described by a simple a-disk model. 

Although we have used a nominal accretion efficiency of 77 = 0.1 and a stress coefficient 
a = 0.1 throughout this paper, these parameters are also uncertain, although by no more 
than a factor of a few. They influence our predictions in only two ways: the total energy 
deposited in the disk and also the radiative cooling time are proportional to S oc a~'^^^r]~^^^ , 
so the total luminosity is independent of either parameter. The inflow time — and thus the 
lifetime of the afterglow phase — is proportional to a~^, so uncertainty in a of order unity 
could vary the total number of observable sources by a similar factor. As discussed above 
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at the end of Section [3l the uncertainty in rh may not be too great, since strong selection 
effects hmit the hkehhood of observing either very small or very large rh. However, if the 
loss cone of scattering stars in the galactic center can be continuously replenished, it may 
be possible to shrink the BH binary sufficiently through three-body interactions alone to 
trigger a GW merger without a circumbinary disk. If so, there may be additional black hole 
mergers without enough surrounding gas to be as bright as the ones described here. At the 
other extreme, there may always be such a large amount of gas driving the merger that no 
appreciable gap is ever formed in the disk, making the afterglow very difficult to distinguish 
from a normal AGN. 

The simple model we use for merger rates is also uncertain by at least a factor of a few, 
but this should only affect the overall number counts of sources, and not directly change 
their appearance. However, if the BH masses in any given merger are not really selected 
at random from the overall mass distribution function (e.g. comparable mass mergers are 
favored due to the shorter time from galaxy merger to BH merger), the number of afterglows 
will increase, since systems with large q tend to fill in their gaps more slowly. Finally, the 
distribution of BH spin orientations is still quite uncertain, and there may be significant 
evolution effects that favor aligned spins with relatively small kicks, in turn reducing the 
average luminosity of the afterglow population. 

Of course, wherever there are theoretical uncertainties, there are great observational 
opportunities. If numerous IR afterglows are successfully detected and positively identified 
as open-gap accretion disks, we can begin to use them as observational tools to probe SMBH 
binary systems. Some of the important questions that may be answered are: what is the 
distribution of BH spins and orientations? what is the range of astrophysically relevant kick 
velocities? what are the merger rates, total masses, and expected mass ratios for SMBH 
binaries at various redshifts? how far out do accretion disks extend and what are their 
properties at large radii? 

We would like to thank Tamara Bogdanovic, Cole Miller, Tim Heckman, and Marta Volonteri 
for helpful discussions and comments. We would also like to thank the anonymous referee for 
their close reading of the text and very constructive comments. This work was supported by 
the Chandra Postdoctoral Fellowship Program (JDS) and NSF grant AST-0507455 (JHK). 
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Fig. 1. — Light curves from a few different characteristic mergers. In all cases, M = 10^ Mq, 
7] = 0.1, and a = 0.1. Solid line: q = 1, l^ick = 1000 km/s, rh = 0.1; Dashed line: q = 1, 
Hick = 300 km/s, rh = 0.1; Dot-dashed line: q = I, Vkick = 1000 km/s, rh = 0.2; Dotted 
line: q = 0.1, Vkick = 300 km/s, rh = 0.1. The step at late times is due to the gap in the disk 
filling in (this happens earlier for smaller mass ratios and higher accretion rates). 
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Fig. 2. — Spectra of perturbed disk at a series of different times for M = 10^ Mq, m = 0.1, 
g = 1, and Vkick = 1000 km/s. The disk brightens at early time, then reddens while dimming 
slightly before brightening again as the inner disk lights up. 
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Fig. 3. — (top) Merger rate of SMBHs in observable universe per year per unit redshift, 
assuming each BH undergoes an average of three major mergers between z = 6 and z = 0. 
(bottom) Total number of potentially observable afterglows per unit redshift. An afterglow- 
ing system is one whose luminosity is dominated by the thermal relaxation of the perturbed 
disk and where the central gap has not yet closed. 
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Fig. 4. — Luminosity distribution function (number per comoving volume per log luminosity) 
for systems in the afterglow phase at different redshifts. As the number of low-mass BHs 
(M < IO^Mq) increases at smaller redshifts, the distribution function shifts to slightly lower 
luminosities. 
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Fig. 5. — Rest-frame spectra of a sample of afterglow systems at any one time (solid curves). 
For reference, the dashed curve is the spectrum of a similar merger system shortly after the 
central gap has closed. 
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Fig. 6. — Contour plot of source counts A^(> S) in the entire sky as a function of observed 
wavelength (vertical axis) and flux (horizontal axis). Also shown are the flux limits for the 
SWIRE (blue diamonds) and COSMOS (red triangles) surveys. More massive systems tend 
to be brighter and last longer, so the A^(> S) relation is flatter than the standard S*"^/^ 
power-law for uniform source distributions. 
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Fig. 7. — Source distribution contour plot as a function of redshift and IR flux, observed at 
(top to bottom) 3.6, 8, and 24^m. The contours are measures of A^(> S) per unit redshift. 



